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I. INTRODUCTION 



The emphasis of this thesis is two-fold. The first point 
is to. gain insight into the characteristics of a spectral 
primitive equation model. The second is to attempt to take 
advantage of the spectral expansions in an initialization 
techn i que . 

A spectral model uses a Galerkin formulation with the 
spherical harmonics as its choice of basis functions. The 
basis functions are eigenfunctions of the Laplacian operator 
in spherical geometry and, thus, are used to advantage for 
global integrations. An immediate consequence of the Galerkin 
formulation is that spatial derivatives of a single wave are 
computed without the usual truncation error present in finite 
difference methods. This means that the spectral model 
contains no linear phase speed errors. 

The non-linear terms are computed via a transform grid 
in physical space (Orszag, 1971). This procedure is respon- 
sible for making the efficiency of the spectral model com- 
parable with finite difference methods. In addition, the 
non-linear terms may be computed without aliasing. This pre- 
vents non-linear instability and conserves invariant integral 
properties. Physical processes are also computed on the 
transform grid. 

The second objective of this thesis is to test a divergent 
initialization for the spectral model. Most static and 
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variational initialization techniques of global primitive 
equation models use a constraint to balance the rotational 
part of the wind with the mass field while setting the 
divergent part to zero. Phillips (I960) has shown that a 
completely non-d i vergent initialization cannot eliminate the 
gravity waves which would be generated and his work suggests 
that the q ua s i -geost rop h i c divergence could tend to suppress 
these gravity wave modes for mid-latitude synoptic scales. 

In addition, if the quas i -geostroph i c divergence had internal 
modes close to those which would be developed by the model, 
such things as the large scale precipitation would be more 
realistic during the first few hours of the model forecast. 

It is an aim of this thesis to develop a diagnostic 
velocity divergence to insert into the fields of an otherwise 
completely balanced system and to test the subsequent effect 
on the spurious gravity wave noise and the developing diver- 
gence. 

A spectral formulation of the primitive equations can 
lend itself to techniques which initialize according to scale 
For example, the q uas i -geost rop h i c divergence might be com- 
puted only for the synoptic scales of the model. Furthermore 
the model conveniently supplies diagnostic information about 
the various scales. 

The final goal of this thesis is to develop a semi-impli- 
cit time differencing scheme and test its effects during the 
first few model hours after initialization. 
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I I . MODEL DESCRIPTION 



The model described here has been developed by Rosmond 
(1977) and the details are presented here for reference 
purposes. Similar formulations have been done by Hoskins 
and Simmons (1975) and by Bourke (1974). The equations for 
an inviscid, adiabatic and hydrostatic atmosphere may be 
written as 

1.1) = -V»(C+f )V - k • Vx(RTVq + a 1^) 

1 . 2 ) = k»Vx(C+f )V - V* (RTVq + - 5 ^) ~ V ((J) + ^ ) 

1 . 3 ) 11 = .7 . ve - 

1.4, |a=-D. 7 . 7 q-|| 

1.5) a||=-RT 

where 

C - vorticity (C = k • V x ^) 

D - velocity divergence (D = V • V) 

T - temperature 

6 - potential temperature 

IT - surface pressure 

^ - horizontal velocity vector 

4 > - geopotential height 

R - gas constant 

Cp- specific heat 
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f - Coriolis parameter 
a - vertical coordinate (a = p/u) 
a - vertical velocity (a = da/dt) 
q - JlniT 

The continuity equation (1.4) may be refined by integrat- 
ing in the vertical and imposing boundary conditions on a; 

• • 

0(0) = 0(1) = 0. Introducing the notation G = V*Vq, we may 
rewr ite Eq. (1.4) as 



by substituting Eq. (1.6) into Eq. (1.4) and integrating in 
the vertical to obtain 



1-6) = -(6 t G) . 



where the over bar denotes a vertical average. 



The vertical velocity, o, may be obtained diagnostically 



1.7) 




o 

By defining a horizontally mean temperature. 



T = T (o) t T'(o, X, 6, t) 



and using the following operator 
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We can now expand the vector quantities in our 
equations in spherical coordinates (X,n) to arrive 



1 .8) 


ae _ 

9t 
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B) 




1.9) 


9D 

9t 


a(B, 


2 2 

-A) -V (r E 


+ (|) + 


.10) 


90 

9t 


-a(ue 


,V6) + 0D - 


. 30 
9a 


.11) 


It 


-(D + 


G) 




. 12) 


34> 


-RT 







where 



A 

B 

G 

E 

U 

V 

X 

n 

y 



, , . , , RT ' Q 

(C+f )U ^ ~ ® 

r 



= (?+f)V - a 



U 9£mT 



= 

2( i V) 

= u cos T)/r 



da 

Ml 

da 
+ V 



aX-nTT 

2 



r 

aJlniT 

9y 



9X 



= V cos n/r 



= longitude 



= latitude 



= sin q 



basic 

at 
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Equations (1.8) - (1.12) are the basic equations used in 
the model. The equations are represented spectrally in the 
horizontal and finite differenced in the vertical. The de- 
pendent variables are written in terms of a triangularly 
truncated series of spherical harmonics: 

X = z xJI y'J 



where the summation is a double sum over m,n for |m| M and 
|m| _< 2, _< M. The separation is such that the coefficients, 
X^, are functions of time and vertical coordinate and the 
spherical harmonics, Y^, are horizontal functions of space. 
The normalization and orthogonality properties of the Y^ 
allow the coefficients to be obtained as follows: 



2tt 1 1 

"" = "X- xjx = 4¥ / f X Y™ d,dX . 

O - I 

The model used for this study has five vertical layers 
staggered as in Fig. I. The even levels carry a and the odd 
levels carry all the dependent variables plus the diagnostic 
variables, U,V. 

The finite difference scheme for the vertical advection 
terms is 



3X 

9a 



2(A0, ) '-*Xk+| 






<Xk - 



This formulation conserves kinetic energy. 
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VERTICAL SCHEME 



gure I. Staggered 5-!ayer sigma coordinate system. 
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Following Arakawa and Lamb (1976), the hydrostatic equa- 
tion (1.12) can be finite differenced as 

“^k+l " ‘f’k ^ ^ ^“k "'^k+l ^k '''k^ 

where 

“k ' ' 6k ■ ‘°k+i''“k>'^''‘^P - ' • 

A boundary condition is formed by integrating the hydrostatic 
equation from o = 0 to I which gives 




where '*'1^® -terrain geopotential. This boundary con- 

dition is then finite differenced as 

1.14) E (l)^(Aa^) - = R Z T^(Aa^) 

k 

Equations (1.13) and (1.14) are then combined to give the 
matrix equation 

1.15) = 1IC]T + 



which is the finite difference form of the hydrostatic equa- 
tion where the variables are now column vectors. 
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The pressure tendency equation (l.ll) can be written as 



la 

3t 



Z + D(^) (AOj^) 

k 



or i n matr ix form as 



1 .16) 



la 

3t 



= -CN] (G+D) 



Similarly equation (1.7) can be written in matrix form as 



1.17) 



b = [Z] (G+D) 



The thermodynamic equation (I.IO) is differenced consis- 
tent with the method employed' in the hydrostatic equation. 
The potential temperature is written in terms of temperature 
and the last term becomes 
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3a 



P. .^k .R/C 

Lcr,.^, (T.^, (- ) p 



2(Aa, ) •- k+l k+l - 
k a 



k+l 



»k'VTk-l‘s— ’ 

Vi 



f^'^k - - 

+ ^ (G^ - G - D). 



For the purposes of the sem i - I mp I i c i t formulation de- 
scribed later in the paper, we separate the temperature into 

* 

its horizontal average, T , and perturbation, T*. The 
matrix, CyD, is defined such that 



* .^k .R/C * 

“k+l 

2(AOk> 



“k" 



. 3|, R/c , 

'k - ’ 'k-l 

^k-l 



2(Aa^) 



a^ = CY]a . 
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The vertfcal ly averaged divergence can be written D = 



Cn!]D, so that 



* 

_[Y]a - ^ D = 



-[Y]a - ^ [N]D 
P 



Substituting for a from equation (1.17) we obtain 
* * 

-CY]a - ^ CN]D = -CY][Z]G -([YXZ] + ^ [N])D 
P P 



= -CS]G + [Q]D . 



Finally the entire term can be written as 



^ la ^ -CsHg+CqDd- 2(Aa, ) ^*k+l ^^k+l ^ ~ "'^k^ 



a. R/C 

t aJT,'-(-!^) 
k k 

Vi 



'^k+l 

, RT, RT' 

k . p . k,^ k ;; 



P P 



The nonlinear terms are computed using the transform 
method suggested by Orszag (1971). The longitudinal direc- 
tion is done using Fast Fourier Transform and the latitudinal 
direction is done with Gaussian Quadrature. The number of 
latitudes, N, and longitudes, M, satisfy: 

N ^ (3J + I )/2 and M >_ 3J + I . 

The number of points are picked to ensure non-alfased results 
for quadratic terms. Certain terms containing a have triple 
products and are not computed alias free. Experiments by 
the author and a previous study (Hoskins and Simmons, 1975) 
indicate that this source of error is negligible for the 
length of the integrations performed in this study. 
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SEMI-IMPLICIT SCHEME 



This section gives details of the sem i - i mp 1 i c i t time 
differencing scheme used in some of the experiments of this 
study. Using the matrices defined in the model description 

section, equations (I.9)-(1.I2) can be written 

2.1) - 5 ^ + + RT q) = 

2.2) - CQ]D = 

2.3) + [nDD = -G 

2.4) <p' = CC]T 



where 



(f) ^ — (j) — (|) 



sfc * 



Fp and F-j. represent the remaining terms in each equation 
which are not explicitly separated out. Following Robert 
et al. (1972), all the terms on the left hand side of equa- 
tions (2.l)-(2.4) are time averaged. We define the following 
averaging operator and represent the local tendencies by a 
Leapfrog finite difference scheme: 



= 



( ) 



t+At 



+ ( 
Y~ 



A-Af 





( - ( 



,t-At 



2(At) 
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Equations (2.1)-(2.4) can then be written as 



2.5) 


6^0 


2 — T-t * — t 

= F^ - V (^^ + RT q^) 




2.6) 


^t^ 


= F^ + CQ] d'^ 




2.7) 




= - G - CnD 




2.8) 




II 

1 1 

o 

LJ 




El i m 1 n a t i n g in equation (2.5) by 


subst i tut i on 


8 ) , we 


arrive at a time averaged set 


of equations: 


2.9) 


d'*’ 


= t(At)(Fj^ - V^([C]T^ + RT* 


q"^)) 


2.10) 




= t(At)(F^ t [Q] D^) 




2.11) 




= q'*’”^''’ -(At)(G t [N] d"^) . 





Substituting equations (2.11) and (2.10) into (2.9) and 
solving for d"*” we obtain a Helmholtz equation for the time 
averaged divergence 

2.12) ([l>(At)^[B]V^)D'^ = + (At)[Fp-V^([C]T''’”^''’ 

+ (At)Cc]F^+RT*q''’"^''’-(At) RT* G] 

* 

where [I] is the identify matrix and CbJ = CCDLOD - RT CN]. 

For a spectral formulation, the equation is trivial to 

solve since the spherical harmonics are eigenfunctions of 

2 — i" 

the Laplacian V . Knowing D , we then compute the prognostic 
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• ..I T-t+At t+At , .. ,, r^t+At 

variables T , q from equations (2.6) - (2.7). D 

is computed from the definition of the time averaging 

op era tor : 



= 2D+ - 



and the vorticity, never entered the sem i - i mp I i c i t 

formalism and is computed from a standard Leapfrog scheme. 
Calculations with the sem I - i mp 1 i c i t model required 5$ more 
computer time per time step. Stable integrations were 
possible with a time step of 60 mins. 
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I V. QUAS I -GEOSTROPH I C DIVERGENCE 



Following the analysis of Simmons and Hoskins (1976), a 
diagnostic divergence Is obtained for the case of quasi- 
geostrophic flow with the full variation of the Coriolis 
parameter over the sphere. The relevant equations become 



3. 1 ) 

3.2) 

3.3) 



§ = - • V(ctf) - <D 

- v^((|)' t RT%) = 0 




VT + § (f 
P 



T 

a 



* 




where 




rotational part of V 
sigma dot at odd levels 
V • (fV(V“^D)) 



Making use of the hydrostatic relationship, equation 
(1.15), we may write equation (3.3) as 



II' = -CC](V^ • VT) + y* b 

where y = [c] (^ - - ) • 

2 9 

Taking V of this equation, Jf of equation (3.1) and 

of eq. (3.2) and substituting we arrive at the following 

A. 

equation for a: 

3.4) = -v‘^C (V^*V(C+f))+ ^D]+CC3(V^^*VT)tRTV || . 
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A second equation for a Is obtained from the continuity 
equation in the form of eq. (I.I7), which when explicitly 
written out, is 



3.5) a 



a(D + G) 




o 

(D + G) do 



o 



where a Is specified on even levels. 

The diagnostic a from equation (3.4) can be used to com- 
pute the divergence using eq. (3.5) provided we specify a 
boundary condition. For the purposes of Initialization, this 
is most conveniently chosen to be 



= - (D t G) = 0 



since this will tend to eliminate external gravity waves. 
Applying this constraint to equations (3.4) and (3.5) we ob- 
tain 

3.6) = -V"^( (V^*V(C+f))t jf^D) + [C]V^'VT 

/ a 

(D t G) da . 
o 

The method of solution will be to compute a using equa- 
tion (3.6) from an initial guess of the divergence. The even 
level a’s are then found by interpolation and a new guess of 
the divergence is obtained by inverting equation (3.7). The 
updated divergence can be substituted back into equation (3.6) 
and the procedure continues iteratively until sufficient 
accuracy is obtained. 
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Equations (3.6) and (3.7) will now be put Into spectral 
form. The equations contain three similar advective terms. 
It will be sufficient to show the spectral form of one of 
them : 

3.8) V^.V(c-^f) =-J(4..C*f) . 



First, the terms are Individually transformed to physical 
space by the summations; 



ic,m 



-.m 



£,m 

dip , , m _m . I mX 

J6,m 

sp'" 

9(?tf ) _ V' pfn n imX ^ 

~5^ 4^ h liT * 

Z,m 



The multiplications in eq . (3.8) are then computed at 

grid points and the result is transformed back to spectral 

space by computing the inner product defined as: 

+ 1 .2tt 



3.9) Jr = x;r 

L 4ir 




J(i(;, C+f)Y^I^ dXdy e <J(i(;, ?+f), . 



The integral is computed using Gaussian quadrature in 
the y direction and Fast Fourier Transforms in the X direc- 
tion. The first term of equation (3.6) can then be 
eva I uated : 

iC(V^'V(^+f)) = V .(fV(V^J)) 
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and expanded in spectral form using the previous computation 
of J^; 

+1 jM • gpM 

M r MM r 7 ^ Q M 

< (Vi|;‘V(C+f)) \> = p'^)pj^dp-fi J (lV)(2 ^>Pl ' 

_ p Sj I 9^ 

where the X integration has already been done. 

Making use of the recurrence relationship 



.pM = M ^ M M 

^£+1 ^£+1 ^£ £-1 



M 



i, 



2 2 
£ - M 

~2 

4£ - I 



we can write the first integrand as 



„MM MMM MMM 



and using the orthogonality of the Legendre polynomials we 

obtain for the first integral 
-+l 



., 0 . / pa -r., <f, 



The second integral can be done similarly using the 
relationship 



8P^ 

, , 2°"^l _ ,0^,. M 0 M JA 

^£ *^£-1 " ^^£+1 *^£+1 
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to obta i n 



+ 1 

>. 11 ) y* ( l-u 



M .pM M M 

2,,^ „M _ ,M ""L+l ,M 

^ „I(ThT 3y L ^ '^L+l L+l ’ \-\ L 



Combining eqs ( 3 . 10 ) and ( 3 . 11 ) we obtain 



3 . 12 ) <#(v^.v(;+f),v”> = ncif!. ^ 



2 

The second order operator X. D can now be formulated. 

-MM 

Using the notation, <*D, Y|^> = Fj^ , we can write by analogy 
with equation ( 3 . 12 ) 



cM _ nrL+l r,M L M ,M 

Fl «C l + |_+| e|_+| J|_+|3 . 



Thus, 



SCO 



I 

l , m 



-M 



Y 



M 

Si 



and 



- ofL+l M pM , L M pM 

<< D, Y^> - </(JfD), Y^_> - £l ^-I ■'l+r^tl ^+1^ * 



M M 

Substituting in for F|__| ^nd we obtain the result 



3 . 13 ) <( 
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Finally, eqs (3.6) and (3.7) may now be fully expanded 



as 






u 

3.15) ^i_ ” “ / 



where 



Rl = Cc: <V^ • VT, Y^> 
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V. 



DESIGN OF THE EXPERIMENTS 



The initial conditions used in the experiments of this 

thesis were developed from a ba roc I i n i ca I I y unstable mean 

flow and a small perturbation in zonal wave no. 6 (m=6). 

The mass field of the mean flow was constrained to be in 

geostrophic thermal balance with the mean wind. The vertica 

profile of the mean wind was linear varying from 4.5 m/sec 

at a = I to 49.5 m/sec at a = 0. The latitudinal structure 
2 

varied as sin (2ri) placing a jet max at n = 45®. The verti- 
cal temperature profile approximated the U. S. Standard 
Atmosphere . 

The perturbation in wave no. 6 corresponded to a maximum 
v-component of .1 m/sec. The model was integrated for 192 
hrs allowing the baro-clinic wave to grow to finite amplitude 
Fig. 2 shows a time series of the vorticity during this 
period. The growth is strikingly linear until 192 hrs where 
there is evidence of the mean flow (m = 0) changing due to 
feedback from the eddies. The I inear growth corresponds to 
a doubling time of roughly 26 hrs for m = 6. The fields at 
192 hrs provided the balanced fields for the initialization 
exper iments. 

Three types of experiments and one control case were 
integrated for 72 hrs of model time. 
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Figure 2. Vorticity time series during the first 192 hrs. 

Vertical axis is logarithmic. Units are sees”'. 
Solid curve is the magnitude of wave no. 0 at 
latitude 60®N. Dotted curve is m=6, short 
dashes are m=2 and long dashes are m=l8. 
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0*261 



The control case consisted of initializing with the 
balanced fields and integrating. The experiments consisted 



of the following: 
Exper i ment A : 

Exper i ment B : 

Experiment C: 



Initialize with the balanced fields but 
no d i vergence 

Initialize with the balanced fields plus 
the qua s i -geost rop h i c divergence 

Identical to Experiment A but with a semi- 
imp I i c i t scheme 



All experiments and the control were first integrated 
with a spectral truncation of wave no. 21. These will be 
referred to as low resolution experiments. High resolution 
experiments were integrated for the control and Experiments 
A and B. These cases were truncated at wave no. 42. Time 
steps of 12 mins for the low resolution and 6 mins for the 
high were used. Experiment C was integrated with time steps 
of 12 and 60 mins. All experiments had a weak Robert time 
filter of 0.1 and no dissipation term was used. 
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VI. RESULTS 



A. CONTROL CASE 

Figures 3-6 show the Initial fields at level 5 for the 
low resolution experiments. Figures 7-10 are the correspond- 
ing high resolution fields. They are characterized by a 
developing frontal zone. The advantage of the high resolu- 
tion is apparent in the temperature gradients across the 
cold front and warm front and in the surface pressure. The 
low in the high resolution surface pressure field is 5 mb 
deeper. The negative and positive areas of the divergence 
and vorticity fields are approximately symmetrical in the 
low resolution fields and quite asymmetric in the high reso- 
lution fields. The region of maximum positive vorticity 
corresponds very well with the maximum temperature gradient. 

The mean fields consist of an indirect cell induced by 
the growing baroclinic eddies. The indirect cell develops 
in order to decrease the vertical wind shear which must remain 
in approximate thermal balance with the decreasing pole to 
equator temperature gradient. The effect is more pronounced 
in the low resolution fields. The high resolution mean 
fields contain a weak indirect circulation and a strong 
d i rect cell. 

The dynamic situation during the period of the control 
case is one of a rapid cascade of energy to smaller scales. 
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Initial high resolution vortlcity for a = .9. Contou 
interval is 5xl0“^sec"l. Mean field is scaled by 10^ 
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The horizontal scale of the frontal zone decreases. Figures 
11-13 show this effect in the high resolution temperature 
fields for times out to 36 hr. Beyond this point, the 
energy builds up at the spectral limit of the model since 
there is no dissipation and the 72 hr fields are non-meteoro- 
logical (Fig. 14 and 15). Because the energy blockage occurs 
at smaller scales, the high resolution case (Fig. 14) is, 
paradoxically, noisier and in greater error. 

B. GRAVITY WAVES 

In anticipation of later sections, we discuss the 
characteristics of inertial gravity waves for the model used 
in this study. The divergence is the most sensitive indica- 
tor of gravity wave motion. As discussed in Hoskins and 
Simmons (1975), equation (2.12) is the sem i - i mp I i c i t analogue 
of the gravity wave equation for a multi~layer model. The 
eigenvalues of the matrix CbD (eq. 2.12) give the permissible 
e i gen ve I oc i t i es of the gravity wave modes. The corresponding 
eigenvectors give the vertical structure. Figure 16 depicts 
the eigenmodes for the experiments in this study. Table I 
gives the periods of the gravity wave modes for zonal wave 
no . 6 in h rs . 



TABLE I 



Externa 1 


Fi rst 


Second 


Th i rd 


Fourth 


3.31 


1 1 . 1 


28.0 


61.7 


1 54 
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Figure 14, 72 hr high resolution temperature for a = .9, Contour 

i nterva I is 5° K, 
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Figure 15. 72 hr low resolution temperature for a = .9, Contour 
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internal divergence modes. Velocities are in m/sec 



The total divergence may be thought of as consisting of 
a sum of the individual eigenmodes: 

4.1) D = Z aj^ 

k 

where the D^. are the eigenvectors. The weighting coefficients, 
a^, may be obtained from the simultaneous equations formed 
by the dot products of the individual eigenvectors with equa- 
tion (4.1). A time series of the coefficients will provide 
useful diagnostic information in later sections. Figures 17 
and 18 are a time series of the external and first internal 
modes for the low resolution control case at a point on the 
60° latitude circle. It is apparent that the coefficients 
contain both the meteorological modes and any excited gravity 
modes. Figure 17 contains no high frequency oscillation and 
Is purely meteorological. Figure 18 contains a meteorologi- 
cal mode with a period of roughly 3 days plus a smaller 
amplitude wave whose period of approximately 1/2 day corre- 
sponds well with the calculated value for the first internal 
mode from Table I. The remaining internal gravity modes 
have periods approaching meteorological values and are diffi- 
cult to d i st i ngu i sh. 

C. EXPERIMENTS A AND B 
I . Low Resol ut i on 

Neglecting the divergence in the initialized fields 
(Experiment A) generated computational gravity waves. Compar- 
ing the pressure tendencies of the control case (Fig. 19c) 
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TIME 



Figure 17. Divergence time series of the external 
mode for the low resolution control 
case. Vertical axis units are sec”^. 
Time is in hrs. 
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L 



with those of Experiment A in Fig. 19a, it is apparent that 
high frequency oscillations of roughly 3 hr period are 
excited. The initial magnitudes of the oscillations are 
approximately 1-2 mb/hr. The oscillations dampen out with 
time as the geostrophic adjustment takes place although the 
dispersal of the gravity waves during the adjustment process 
is hampered by the longitudinal periodicity and the symmetry 
of the fields about the equator. 

Figure 19b shows the pressure tendencies with the 
q ua s i -geostrop h i c divergence Inserted into the initializa- 
tion fields (Experiment B). The highest frequency oscilla- 
tions have been reduced over Experiment A by a factor of 
ap p rox i mate I y 2 . 

Comparison of the divergence modes in Figs 20, 21, 22 
allows us to be more quantitative. The figures have the 
control case subtracted out and contain the gravity wave 
noise generated by the approximations of each experiment. 

The noise generated In the internal modes by the non-d i vergent 
initialization is an order of magnitude larger than the ex- 
ternal mode and is generally 25$ of the total divergence. 

The figures show that the external mode of Experiment B is 
about 68$ of Experiment A. The first internal mode shows 
more noise in Experiment B initially and then decreasing to 
50$ and the second internal mode shows a decrease of about 
50$ over Experiment A. In general the effect of the quasi- 
geostrophic divergence is to reduce the gravity wave noise. 
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Figure 22. Gravity wave noise in the second internal mode - otherwise 
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Figure 21, Gravity wave noise in the first internal mode - otherwise 
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Figure 22, Gravity wave noise in the second internal mode - otherwise 



2 . 



High Resolution 



Figure 23 shows the q u a s i -geost rop h i c divergence at 
level 5 which was used for Experiment B. A qualitative com- 
parison with the actual model divergence (Fig. 10) shows the 
phases to be in agreement but that the q ua s i -geost rop h i c 
approximation is an overestimate for the perturbation fields. 
The zonal mean fields do not compare well. The quasi- 
geostrophic divergence shows a strong indirect circulation 
which is displaced poleward relative to the model divergence. 
There is no evidence of a direct cell. 

Figures 24-26 depict the noise generated in the 
divergence modes for Experiments A and B. Generally they 
show that the gravity wave noise of Experiment B is 66 %~ 12 % 
of that in Experiment A. The effect of the q ua s I -geo s t rop h i c 
divergence is to reduce the noise. However, this reduction 
is less than In the low resolution experiments. 

In view of the fact that the high resolution experi- 
ments were less successful than the low resolution and that 
quas i -geostrop h i c formulations are valid only for small 
values of the Rossby number, it is interesting to compare the 
initial model divergence of the control case with the quasi- 
geostrophic divergence in a spectral manner. Referring to 
Fig. 27, it is seen that there is good agreement between the 
model and q ua s i -geos t rop h i c divergences for zonal waves 
corresponding to m = 6, 12, 18. The overestimation by the 

quas i -geostroph i c divergence noted earlier is primarily in 
wave no. 6. The mean fields (m=0) and waves corresponding 
to m = 24,30,36,42 do not compare well. 
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Figure 23. Quas i -geostroph I c divergence (high resolution) for 
Contour intervai is IxlQ-^sec-l. Mean fieid is sc 
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Figure 24. Gravity wave noise in the external mode for high resolution 
Experiment A (24a) and Experiment B (24b). Vertical scale 
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Figure 25. Gravity wave noise in the first interval mode - otherwise 
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Figure 26 . Gravity wave noise in the second interval - otherwise 
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Figure 27. Initial divergence coefficients for the control case 

(dotted line) and Experiment B (dashed line). Verticai 
axis is latitudinal wave number (lat. index of the 
spherical harmonics). Horizontal scale is sec"^. 



One motive for initializing with the quas i -geost ro- 
phic divergence is that, presumably, the "true” divergence 
would be approximated at the start and that the model would 



thus develop this divergence field earlier than a non- 
divergent initialization. Figures 28 and 29 address this 
point. They show that at 6 hours into the integrations, zonal 
wave numbers 6, 12, 18 do, in fact, approach the control 

divergence slightly faster. There is almost no difference 
for the higher wave numbers and the mean field shows poor 
adjustment in either case. This might be expected for the 
mean field since the geostrophic adjustment process is very 
inefficient at this scale in a bounded region. On the other 
hand, the higher wave numbers adjust very rapidly and after 
6 hrs, they are in good agreement with the control case. 

In light of the results showing good agreement in 
the initial quas i -geostroph i c divergence with the control 
divergence for zonal wave number 6, 12, 18 and bad agreement 

for other scales, it is natural to modify Experiment B such 
that only scales corresponding to m = 6, 12, 18 are used 

in the qua s i -geost rop h i c divergence. Figure 30 shows the 
time series of the external divergence mode for this experi- 
ment. Figure 31 shows the spectral adjustment of the diver- 
gence at 6 hours. The corresponding results for the unmodi- 
fied Experiment B are Figs. 24 and 29. The results are 
almost identical to the unmodified q ua s i -geos t rop h i c diver- 
gence results. 
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Figure 28, 6 hr divergence coefficients for the control case (dotted 

line) and Experiment A (dashed line) - otherwise as in 
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Figure 29. 6 hr divergence coefficients for the control case (dotted 

line) and Experiment B (dashed line) - otherwise as in 




TIME 



Figure 30. Time series of the gravity wave noise 

in the modified Experiment B. Vertical 
scale is sec-l. Time is in hrs. 
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Figure 31. 6 hr divergence coefficients for the control case (dotted 

iine) and modified Experiment B (dashed line). Vertical 
axis latitudinal wave no. Horizontal scale is see"'. 



D. EXPERIMENT C 



As discussed in Hoskins and Simmons (1975), a linear sta- 
bility analysis of the sem i - i mp I i c i t scheme shows that the 
time step is limited only by the Rossby waves. The effect 
of the scheme is to slow down the faster moving gravity wave 
modes so that the integrations do not amplify in time. 

Slower moving atmospheric motions are less affected. Table I 
presents sample calculations of the distortion of the periods 
of the gravity waves present in these experiments for two 
time steps. These may be compared with Table I. The general 
behavior is that the periods increase with increasing time 
steps. 



TABLE I I 



At 


Externa 1 


First 


Second 


Th i rd 


Fourth 


12 mins 


3.31 


11.14 


28.0 


6 1.7 


154.0 


60 mins 


5.7 


12.2 


28.5 


6 1.9 ' 


154.1 



Two points will be investigated in this experiment. The 
first is that the sem i - i mp I i c i t scheme should tend to smooth 
out high frequency oscillations, particularly with a large 
time step. The second is that the sem i - i mp I i c i t method 
might possibly hamper the geostrophic adjustment process 
since it slows down the gravity waves which are the mechanism 
for the adjustment. 
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Figure 32 shows the pressure tendencies using the semi- 
implicit method with time steps of 12 and 60 mins. Compar- 
ing this with the corresponding results for the explicit 
integrations (Fig. 19a), it is evident that the effect of 
the sem i - i mp I i c i t method is to smooth out the 3 hr oscilla- 
tions. Figure 33 shows the relative effect of a small vs. 
large time step in the sem i - i mp I i c i t method on the external 
gravity waves. The larger period of the 60 min time step is 
evident. The larger time step initially amplifies the com- 
putational noise, but after 24 hrs it is considerably damped 
The smoothing In the pressure tendencies is postulated to 
come about from the decreased velocity of the external 
gravity mode assuming an advective time scale. 

The second point of this experiment is shown in Figs. 34 
and 35. Comparison of the divergence coefficients for time 
steps of 12 and 60 mins shows no difference at 3 hrs. The 
adjustment of the divergence to the control case is apparent 
ly the same for either time step. 
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Figure 32. 3q/3t for Experiment C with a time step of 12 mins (32a) 

and 60 mins (32b). Vertical scale is sec“l. 
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Figure 33, Gravity wave noise in the externai mode for Experiment 
with time steps of 12 mins (33a) and 60 mins (33b). 




mPLITUDE 





Figure 34 . 



3 hr divergence coefficients for the low 
resolution control case (dotted line) and 
Experiment C with a 12 min time step 
(dashed line). Vertical axis Is latitudinal 
wave no. Horizontal scale is sec"! 
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Figure 35. 3 hr divergence coefficients with a 60 

min time step - otherwise as in Fig. 34. 
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SUMMARY AND CONCLUSIONS 



VII. 

The results of this study are based upon the divergence 
developed by a baroclinic wave on a uniform earth. Initiali- 
zation without the divergent wind component generated spurious 
inertial gravity waves. The internal state of the initial 
fields had nearly zero mean divergence. Thus, the external 
mode of the generated gravity waves was an order of magni- 
tude less than the total divergence and was not significant. 
However, the internal gravity wave modes generated were signi- 
ficant. An important source for the external mode, namely 
topography, was neglected. 

Initializing with the q ua s i -geost rop h i c divergence 
generally decreased the spurious gravity wave noise and aided 
the ensuing integrations to develop a divergence closer to 
the control case, but both effects were small. Furthermore, 
the advantage of the q ua s i -geost rop h I c divergence decreased 
when smaller scales were included in the integrations. In- 
serting the qua s i -geostroph i c divergence into only synoptic 
scales gave no advantage. However, the effect of this experi- 
ment may have been reduced by the fact that there were no 
planetary scale waves. 

Using a semi-impi icit scheme with the same time step as 
an explicit integration had no effect on the divergence modes 
but did smooth out oscillations in the pressure tendencies. 
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A larger time step Initially increased the amplitude and 
period of the fast external gravity waves and had less effect 
on slower moving motions while still stabilizing pressure 
tendencies. 

In general, the effect of a divergent initialization for 
a global primitive equation model is small. This same con- 
clusion has been reached by two similar studies - Houghton 
et al (1971) and Dey et al (1975). This study, however, has 
neglected any physical processes such as precipitation as 
well as orography and the effect of an initial divergence In 
these cases may well be more important. 
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